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Abstract 

We consider a family of time dependent dephasing Lindblad gener- 
ators which model the monitoring of the instantaneous Hamiltonian of 
a system by a Markovian bath. In this family the time dependence of 
the dephasing operators is (essentially) governed by the instantaneous 
Hamiltonian. The evolution in the adiabatic limit admits a geometric 
interpretation and can be solved by quadrature. As an application 
we derive an analog of the Landau-Zener tunneling formula for this 
family. 

Lindblad generators describe the quantum evolutions of finite systems 
coupled to a memoryless bath [I]. They give a useful phenomenological 
description of thermalization and decoherence (2J [3] . We consider a family of 
time-dependent dephasing Lindblad generators, first introduced in [I] , which 
models the continuous monitoring of the instantaneous energy. This is the 
case, for example, in the Zeno effect [5]. The family can be defined for 
arbitrary dephasing rate, however, its physical interpretation in the strong 
dephasing regime requires some care [6]. 

The family of dephasing Lindbladians that we consider is denned in such 
a way that any instantaneous stationary states of the Hamiltonian are also 
instantaneous stationary states for the Lindbladian. This makes the family 
special and non-generic. (Generic Lindbladians have a unique stationary 
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state). In the adiabatic limit, the evolution generated by this family can be 
solved by quadrature and admits a geometric interpretation. The qualitative 
features of the dynamics differ from the corresponding dynamics of generic 
Lindbladians [7] reflecting the special character of the family. 

In 1932 Landau [8] and independently Zener [9] and Majorana [10] found 
an explicit formula for the tunneling in a generic near-crossing undergoing 
unitary, adiabatic evolution. Here we describe the corresponding result for 
the non-unitary case associated with dephasing Lindbladian. The solution 
appears to be the simplest generalization of the Landau Zener problem which 
is still exactly soluble. 

The influence of dissipation and noise on the tunneling of a two level 
system has been extensively studied in the physics literature [TTJ [J2J, [T3"| IT3]. 
The results presented here differ in one or both of the following aspects: 
First, in the framework: We assume a Lindbladian, and do not attempt to 
derive an effective dynamics from a (unitary) model of a bath or a (unitary) 
model of stochastic noise source. Second: The Lindbladians are, as stated, 
of the dephasing type. Our results contribute to the mathematical physics of 
Lindblad operators. This approach has the virtue that the adiabatic dynam- 
ics can be solved by quadrature and does not rely on the assumption of weak 
dephasing which one normally needs to make when modelling decoherence 
with a unitary bath or noise. 

In the limit of weak dephasing, the tunneling formula we derive can be 
compared with results of [13] for Zener tunneling due to a dephasing noise. 
The two formulas have the same functional form up to an overall constant 
which is left undetermined in [13]. 

Since tunneling is dominated by the near crossing dynamics, the universal 
aspects of near crossing in a two level system are captured by a Hamiltonian 
that depends linearly on time. By an appropriate choice of basis and of the 
zero of energy the relevant dynamics is governed by the Hamiltonian 



where e > is the adiabatic parameter and go > is the minimal gap. 
The tunneling probability T is the probability of a state, which originates 
asymptotically on one eigenvalue branch, to end up in the other at late times. 
The formula Landau and Zener found [15] for this Hamiltonian is: 




(1) 
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The singularity of the limit he — > reflects the singularity of the adiabatic 
and semiclassical limits, and their coincidence in this case. 

The universal aspects of tunneling for near crossing in an open system 
described by a dephasing Lindbladian is described as follows. The adiabatic 
evolution of the density matrix p is governed by 

hep = C s (p), (e > 0) (3) 

where the slowly varying parameter s = et, having the physical dimension of 
an energy, is viewed as the slow clock. C s is the changing Lindblad operator 

C(p) = -i[H,p] - h 1 (P_ P P + + P +P P^ (4) 

H is the Hamiltonian, which for a generic near crossing is given in Eq. (PQ); 
P± = |±) (±| are the two spectral projections of H; finally, 7 > is the 
dephasing rate [16]. 7 = is the case considered by Landau and Zener. 
In both cases transitions between the ground and the excited states only 
occur because the generator of the dynamics depends on s. The tunneling 
probability 

T = tr(pP + )(oo), (p(-oo) = P_(-oo)) (5) 

is the error in fidelity of the ground state. 

The adiabatic tunneling formula with dephasing, which we shall derive 
below, is [T7] 

where Q is the algebraic function (shown in the figure) 

Q(x) = 1 ^ + v / TT^ ^_ 

Few remarks about this result are in order: 



T = ^ L 2 Q[-)+0(e 2 ), (6) 



The adiabatic limit means that v he is the smallest energy scale in the 
problem and in particular, e /ry 2 . When this fails, the error terms 
in Eq. need not be small compared to the leading term. 



When the dephasing is weak, fry -C go, Eq. (JS)) reduces to 

37r e^h 2 
16 



(8) 
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Figure 1: The function Q(x). The argument is the ratio of dephasing rate to 
the minimal gap. The function has a maximum at x = 1.13693 



This term has the same form as one of the tunneling terms found by 
Shimshoni and Stern [13] in a (different) model where the two-level 
system is dephased by noise. The method they use cannot give the 
overall constant 37r/16 [IS] , nor does it allow investigating the full range 
of hj/g . 

• When hq ^> g Eq. dSJ) reduces to 



This may be understood as a manifestation of the quantum Zeno effect 
[5]: The dephasing term in the Lindblad generator can be interpreted as 
monitoring the state of the system at rate 7. This suppresses transitions 
between the states. 

Eq. follows from, and is a special case of, a more general and basic 
formula for the tunneling when the adiabatic evolution takes place on a finite 
interval of (slow) time [s , Si] and one also allows 7(5) to be time-dependent 




ds + 0(e 2 ), 



(10) 



Here g(s) is the instantaneous gap in H(s), 

g\s) = s 2 + g 2 . 



(11) 



The positivity of the integrand in Eq. (Till]) when 7 > makes the tunneling 
irreversible. This changes the characteristics of the e dependence of T from 
exponentially small in Eq. (J2]) to linear in Eq. ((6]). The Landau-Zener formula, 
Eq. ([2]) , is buried in the error terms of Eq. ffTO]) . 

Eq. f flUj) reduces the tunneling problem to integration. In the case where 
7 is constant and s runs from —00 to 00 the numerator Eq. ffTUl) is simply 

tr(P + P 2 P + ) = JL- y (12) 

Elementary algebra then leads to Eq. ([6]) with 

/oo 
(t 2 + l)" 2 (t 2 + l + x 2 )- 1 dt. (13) 
-00 

The integral can be evaluated explicitly to give Eq. (J7J). 

The key idea behind the derivation of the adiabatic tunneling formula, 
Eq. ( FlOl) . is a geometric view of the spectral projection as adiabatic invari- 
ants. The evolution of observables is governed by the adjoint of the Lindblad 
generator, £*, (this is the Heisenberg picture), where the adjoint refers to the 
Hilbert-Schmidt inner product. In particular, the adjoint of the dephasing 
Lindblad operator of Eq. (j3J) acting on the observable A is given by (from 
now on we set h = 1) 

C*(A) = i[H, A] - 7 (P_AP+ + P+AP_), (14) 

It differs from Eq. (J4]) by the replacement of % by — %. As we shall now see an 
instantaneously stationary observable A(s) G Ker(£*) that has no motion in 
Ker(£ s ) is an adiabatic invariant. More precisely, 

Theorem 1. Let A(s) be an observable which lies in the instantaneous kernel 
of C* s , i.e. 

C* s {A( S ))=0 (15) 
and suppose that, in addition, the linear equation 

A(s)=C:(X(s)) (16) 

admits a solution X(s). Then one has 

tr(A( S )p £ ( S ))\Z = etT(X(s)p £ (s))\l - e H tr(X(s)p £ (s)) da, (17) 

J SQ 
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where p e {s) is a solution of the adiabatic Lindblad evolution. A(s) is an 
adiabatic invariant in the sense that its expectation is conserved up to a small 
error, 0(e), given by the right hand side of Eq. ( fT/j ) whereas the change in 
a generic observable is 0(e~ l ) and in the Lindblad generator is 0(1). 



The identity, Eq. ffTTjl . readily follows from 

^ tr(A(s)p £ (s)) = tr(A(s)p £ (s)) + tr(A(s)p £ (s)) 

= tr(C:(X(s))p £ (s)) + e- 1 tr (A(s)C s (p £ (s))) 
= tr(X(s)C s (p £ (s))) + e~ l tv(jO* s (A(s))p £ (s)) 
= etr (X(s)p £ (s)) (18) 

and integration by parts. 

Eq. (ITEj) may be interpreted as a condition that A(s) undergoes parallel 
transport: The equation has a solution provided A(s) G Range (£*) which is 
the case if A(s) has no motion in Ker (C s ). 

It is straightforward to verify that the instantaneous spectral projections 
Pj(s) of a dephasing Lindblad generator are adiabatic invariants in the sense 
of the theorem. Evidently, £*(P + (s)) = 0. Moreover, Eq. ffTBT) is solved by 

X(s) = -iY PkP+Pl (19) 

with e± the two eigenvalues of H. To see this note first that X(s) is purely 
off-diagonal [19] by construction and so is P+, namely 

P+ = P-P+P+ + P+P+P- (20) 

This follows form P + = P+, which implies P + = P+P+ + P+P+ and in turn 
P±P + P± = 0. The equality of the off-diagonal components of Eq. f|T6l) follows 
from 

C*(P k APj) = i(e k - ej + %i)P k AP 3 , (k,j = ±,k^j). (21) 

The probability of leaking out of the instantaneous ground state is given 
by Eq. (ITTj) with A(s) = P+(s). Eq. (fTOjl then follows by appealing to the 
adiabatic theorem [7] which allows to replace the instantaneous state by the 
instantaneous projection on the right hand side of Eq. ([Tj 



P_(s) + 0(e), (22) 
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uniformly in Sq, s. The rest is simple algebra. 

For the convenience of the reader we include a proof of Eq. (1221) . Let 
U £ (s,sq) be the propagator for the differential equation ([3]), whence p(s) = 
U £ (s, so)-P-(so). We recall that the solution of its inhomogeneous variant, 
x = + y, is given by the Duhamel formula 

x(s) = U £ (s,s )x(s )+ f U £ (s,s')y(s')ds'. (23) 

J so 

The remainder to be estimated, r(s) = p(s) — P_(s), satisfies r = e~ 1 C s (r) — 
P_, because of £ S (P_) = 0. Before applying Eq. f[2~3"j) . we observe that the 
equation P-(s) = C s (X(s)) admits Eq. ffT9l) as a solution upon replacing i, 
P + by —i, P_. The differential equation thus becomes 

(r - eXy = e- l C s {r - eX) - eX, 

resulting in 

r(s) - eX(s) = U £ (s, s )(r(a„) - eX(s )) - e f U £ (s, 8')X( 8 ')ds'. 

J SO 

Since C s is dissipative, i.e. tr(p£ s (p)) < 0, we obtain ||C/ 6 (s, so)|| < 1, (s > 
s ). Together with r(s ) = we conclude that r(s) = 0(e), as claimed. The 
uniformity follows from decay: X(s) = 0(s~ 3 ), X(s) = 0(s~ 4 ), (s — > ±oo). 

In conclusion: We have introduced a class of adiabatically changing de- 
phasing Lindblad operators which allowed us to calculate the tunneling in 
a generic two-level crossing and extend the Landau-Zener tunneling to de- 
phasing Lindbladians with arbitrary dephasing rate. Dephasing makes the 
tunneling irreversible and so fundamentally different from tunneling in the 
unitary setting. This irreversibility is responsible for the difference in the 
functional form of the tunneling formulas. 
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